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Abstract: We study the growth of perturbations in an expanding Bianchi type-I metric 
that evolves according to an energy density that includes dust and a cosmological constant. 
Assuming an epoch where the cosmological constant is subdominant, we find that, for a 
reasonably large set of initial conditions, the metric fluctuations grow fast enough to make 
the metric inhomogeneous before the cosmological constant becomes the dominant form of 
energy. We have examined values for the cosmological constant that are in the interval 
10 10 GeV to 10 16 GeV. 
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1. Introduction 

Inflation offers so far the best description for the power spectrum of density fluctuations with 
its slight deviation from scale invariance. The claims that inflation explains the "homogeneity 
and flatness" are more questionable. There is an extensive literature on trying to ascertain 
under which conditions for the initial metric and the energy-density will a space develop 
a period of inflation [1]-[13]. Although there is wide spread belief that space-times with a 
cosmological constant isotropize at late times, there is no complete formal proof [14]. One 
result on which there seems to be general agreement is Wald's proof [3] that all expanding 
anisotropic but homogenous models in the presence of a positive cosmological constant, with 
the exception of some Bianchi-IX models, evolve toward a de Sitter solution. In the case 
of the Bianchi-IX this result is also true provided the cosmological constant is sufficiently 
large compared to spatial curvature terms. Recent work has also been directed toward un- 
derstanding the possible experimental signatures of an initial universe that is homogeneous 
but anisotropic [15] - [20]. 

In this work we consider an initially expanding type-I Bianchi metric that evolves accord- 
ing to an energy density not yet dominated by a cosmological constant. This metric is one of 
the special cases studied by Wald [3]. According to his work, if the metric maintains its homo- 
geneity until the cosmological constant becomes the dominant form of energy, it is expected 
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to asymptote to a de Sitter cosmology. The fluctuations of this metric, however, will make 
it evolve from an anisotropic but homogeneous metric into an anisotropic and inhomogenous 
one. The question we seek to answer is: will this growth be fast enough as to invalidate the 
assumption of homogeneity before the cosmological constant becomes the dominant form of 
energy? We will see that this is indeed the case for a reasonably large set of initial conditions, 
when we assume for the cosmological constant to be in the interval between 10 10 GeV and 
10 16 GeV. 

In the next section we will summarize the known results about the time evolution of the 
type-I Bianchi models and its fluctuations, and in section 3 we will present our results. 

2. Evolution of Type-I Bianchi and its fluctuations 

The type-I Bianchi metric is of the form 

ds 2 = -dt 2 + e 2a (e 2 % dx { dx j 

where a = a(t) and characterizes the volume expansion, and fyj 
traceless metric that describes the anisotropy, 

3 

(e 2 % = 5 t ,e 2 ^ J2& = ° ( 2 - 2 ) 

i=l 

The volume factor is only a function of a, y/—g = e 3Q . The equations of motion are described 
in Appendix-A, in this section we will only present the solutions for the dust background, that 
is for an energy momentum tensor whose only non- vanishing component is Too = w o e~ 3a - In 
this model we have made the simplifying assumption that the unperturbed energy-momentum 
tensor is isotropic, placing all the anisotropy on the initial metric configuration. We have cho- 
sen this energy-momentum to take advantage of the existing analytical solution. In addition, 
this choice of the energy-momentum tensor will not exacerbate the anisotropy, as the solu- 
tion below shows. Therefore, if non-trivial results are obtained for this example, we expect 
this behavior to generalize to more anisotropic forms of the unperturbed energy momentum 
tensor. 

Pij = hj u + Cij (2.3) 

is a solution provided we identify 

du = e~ 3a dt (2.4) 

and the matrices bij and c^- are traceless, diagonal and constant. The functions in the metric 
evolve as: 

e^ = ^t{t + t b ) 



(2.1) 

Pij(t) is a 3 x 3 diagonal 
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(2.5) 




In this solution, b 2 = ^hjbij, % = 4b/wo, s« = bu/b and the Cy have been set equal to zero. 
The Si determine the shape of the anisotropy and satisfy the constraints: 
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J> 2 = 6, ]>> = ( 2 - 6 ) 




From these solutions we can observe that at very late times t >> ij, the effect of the 
anisotropy is washed out and the metric evolves as an isotropic FRW cosmology in the pres- 
ence of dust. It is for t « % that the model shows its maximum anisotropic behavior. 
Indeed, although the overall volume factor, \J zz g, increases, this is not necessarily the case 
for the direction-dependent scale factors. Labeling the Sj exponents in increasing order, the 
constraints (2.6) require that 



(A note of caution: in the remaining of the paper we might use a different order for these 
indices, and S3 can be a negative number.) The different components of the metric in the 
limit t << tb, evolve as: 



With the exception of the point (#1,52,53) = (—1,-1,2) the other solutions to (2.6) give 
cosmologies on which, at least, one direction is contracting. 

A complete Bardeen-type analysis of the perturbations to this metric and identification of 
gauge invariant quantities has been recently done by Pereira, Pitrou and Uzan [15]. Because 
our work was started before their work was published, we use an earlier work of Perko, 
Matzner and Shepley (PMS) [21], where the equations are solved in the synchronous gauge 
instead. Since the synchronous gauge leaves a remaining gauge invariance, care will have to 
be taken to assure that the results that we use are not a gauge artifact but actually reflect a 
physical behavior. In our case, since we restrict our study to dust, the synchronous gauge is 
co-moving and the gauge is completely fixed, removing any source for concern on this issue. 
Earlier work by Noh [24] ( see also [26]) who repeated Perko, Matner and Shepley computation 
in the co-moving gauge, showed agreement in the quantities used in this work. 

In the synchronous gauge, the perturbed metric can be written as: 



-2 < si < -1 < s 2 < 1 < s 3 < 2 



(2.7) 



2 = g = e 2«+2ft „ t 2(l+ Si )/3 



(2.8) 




(2.9) 



The energy momentum tensor for a perfect fluid in its rest frame is 



T^ u = V9ixv + ip + ^)u^u, 



(2.10) 
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where u M = 
written as 



(— -^,0), = (— • v /7,0). The perturbed density, pressure and fluid velocity are 



p = p + 5p 
u'v = u fl + 5u» 



(2.11) 



With our gauge choice 5u° = 5uq = 0, while 5u l / 0. 

The details of the perturbations analysis can be found in Appendix B. The metric tensor 
decomposes into two scalar, two vector and two tensor degrees of freedom. Our restriction 
to a matter fluid, although maybe not as realistic as a radiation fluid, greatly simplifies 
the calculations since it is possible to operate in a comoving (in addition to synchronous) 
gauge. Also, the wave vector of the Fourier components of the modes is for simplicity fixed 
along one of the spacial directions (specifically, the third direction, k = (0, 0, ks)). With 
these assumptions, one of the tensor modes is free and the vector modes decouple from other 
perturbations. The equations of interest to our analysis involve a scalar, the density contrast 
(the other scalar modes depend on it), and the remaining tensor degree of freedom, which are 
coupled to each other 



+ 



{^){(FK 2 t 2 ){t + t b ) 2 + \s 2 t 2 } 

(t + h) 



st (t + t b ) 



2 + 2\ 



'-(FK 2 t 2 ) ' 233 



(t + t b ) 2 {FKH 



dS 



T tb + 3 i2t + tb) ) + h dt 



l -{ Sl -s 2 ){FKh 2 ) 



Ml) " + 7" + ^ 



-A 4 



(2.12) 



d271 + ^K d 4 + r ] K 2 F = \ Sl -s 2 ) 



t b d5 



dt 2 ' t(t + t b )dt ' ''" ^ 3^ " ZJ t(t + t b )dt <y2A3 " > 

where S and rj are respectively the Fourier components of the density constrast and of the 
tensor mode (the Fourier index k has been suppressed), 



F = 



1 



2 2£3 

3 /t + t b N 3 



t{t + t b )J V t 

Depending on the value of S3 we can distinguish two cases: 



/ 4 \ 

and K = kp 1 I 



\3ujq J 



(I) S3 = ±2 (i.e. in the axially symmetric case s\ = S2), the density contrast and the tensor 
mode decouple and the solutions are 
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4 + s 



2 (t+t t )3 



for S3 = 2 
for S3 = —2 



(2.14) 



where A and B are constants of integration. 



(II) S3 / ±2, analytic solutions can be found in the long (FK 2 t 2 <C 1) and short (FK 2 t 2 » 
1) wavelength limits. For long wavelengths we have 



A 2(2-s 3 ) / t 

5 = - + B + Ct^ 1 1 + Dln- 
t \ t b 



(2.15) 



and for short wavelengths 



(l-2s s ) 

5 = At 3~~^ + B 



(2.16) 



These solutions are derived in Appendix B. 
3. Study of Perturbations 

In the previous section we derived the time dependence of the matter perturbations. In 
this section, we will introduce a useful parametrization for analyzing their growth. We will 
not be interested in studying the axially symmetric cases S3 = ±2 for two different reasons. 
When S3 = +2 the perturbation doesn't grow in time (2.14) and hence it can never become 
important. For S3 = —2 the growth is wavelength dependent and would require a separate 
investigation. In the limit t <C i& and S3 / ±2, for the perturbation that is outside the 
horizon, the fastest growth is given by: 



while inside the horizon 



2(2-33) 
~ t 3 



(1-2.3) 
~ t 3 



So as long as — 2 < S3 < \ the perturbations will grow both inside and outside the horizon. 
Let's introduce the following parametrization 



5 = 10" 



Sp 



tp 



t b = 10 n n P 



(3.1) 
(3.2) 
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where tp = 5 x 10 _44 s is the Planck time, n and n b are positive integers and s stands for the 
power s> = 2 ( 2 ~ S3 ) outside the horizon or for s< = - | Sa - ) inside the horizon. 
The time t\ at which the perturbation becomes of order one is 

6\ t = tl = 1 => ti = t b 10- n » + ^ (3.3) 

For consistency, we shall impose the condition t\ <C t b since the perturbation analysis 
carried out so far is for times well before the anisotropy-FRW transition. It then follows that 
n and n b are constrained by 

n b > - (3.4) 
s 

M 2 

The total energy density for a pressureless fluid at t <C t b is w ~ . From the preceeding 
equations it follows that 

u\ t=tl = Ul ~ M P lO-~ {nb+ ^ (3.5) 

Let's parametrize the scale of inflation A = MplCT 9 , where q is a positive integer. For the 
matter perturbations to reach order one well before the universe starts to inflate, ui* S> A, 
which translates into the constraint n& + — < 4q. This last equation, combined with (3.4) 
produces the condition - < 2q, which will be later employed in order to derive a possible 
range of values for t\, t b , uj±. 



3.1 Growth of perturbations at different wavelengths 

Although the growth of the scale factor, for t <C h, is different for different directions, 



a t (t) = e a+ ^ 



. q 1 1 + s,- 1 ~ s i 

(^) 3 *-V t«t b 



the horizon scale has the same growth for all of them. 



Hi{t) 



) j t<£t b 



• i+gj 

3 

37 



As already seen in the previous section, and in analogy with the isotropic case, the pertur- 
bation growth depends on the size of the physical wavelength in regards to the horizon scale. 
Indeed, the critical quantity FK 2 t 2 is nothing but 



FK 2 t 2 = constant 



( h 

\a 3 H 
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where H ~ l/f. The product a^H evolves with time as: 

a 3 H(t) ~ ( ^ * « «6 

i t~T t > t 6 

Thus wavelengths that are inside the horizon at early times (close to the Plank scale) will 
remain inside the horizon at later times. The wavelength A decreases as a function of time 
until horizon crossing for — 2 < S3 < — 1 and increases if — 1 < S3 < ^. 

It is important to search for constraints to the parameters involved in the perturba- 
tion analysis. Let's begin from the time t& of transition between anisotropy-dominated and 
FRW-like universe. The universe appears isotropic and homogeneous up to the highest ac- 
cessible redshifts, therefore an upper bound for % would be around 13.7 x 10 9 yr. The lower 
bound would obviously be reprensented by the Plank time. The scale of inflation has been 
constrained from above to 10 16 GeV by the gravitational wave experiments [28]. This limit 
translates into the condition 3 < q. 

No lower observational bound on the initial values of the matter perturbations can be 
determined since the perturbations would grow to order one and the outcome of non linear 
theory would be independent of initial conditions. 

Finally it is important to remember that the expressions derived for the matter pertur- 
bations are valid in two special limits, short and long wavelengths, and the values of q have 
to be in the intervals derived in this section. The lower bound on q and upper bound on % 
can be used to further constrain n and n&, as follows 



q = 3 => n < 6s, - < n b < 12 , (3.6) 

s s 

t b < 10 6 % =>■ n b < 60. (3.7) 

In the Table 1 the parameters have been calculated under the condition that the (inside 
the horizon) matter perturbations reach order one before the time corresponding to the ear- 
liest possible beginning for inflation. 



Table 1: Range of variation of parameters for S3 = — v/3, q=3 



n 


n b 


5 P 


t\ (in units of tp) 


1 

(oji) 4 (in units of Mp) 


% (in units of tp) 


8 


5.3 to 6.6 


10- 8 


2(10 5 ) 


10~ 3 to 2(10~ 3 ) 


h to 4(10 6 ) 


6 


4.0 to 8.0 


10~ 6 


10 4 


10~ 3 to 10~ 2 


h to 10 8 


4 


2.7 to 9.3 


10" 4 


5(10 2 ) 


10- 3 to 4(10- 2 ) 


h to 2(10 9 ) 


2 


1.3 to 10.6 


10~ 2 


20 


10~ 3 to 0.2 


ii to 4(10 10 ) 


10- 1 


0.06 to 11.9 


0.8 


1 


IO- 3 to 0.9 


h to 7(10 n ) 
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As expected from the definition of n (3.1), the perturbations reach order one at a smaller and 
smaller time as n is given smaller and smaller values. Also as n approaches zero, the range 
of allowed values for rib becames larger and larger, and consequently the same happens for 
the range of possible values of tb- If the upper bound on the inflationary energy scale was 
pushed down by a few orders of magnitude, the range of possible values for n and rib would 
be broader and the same conclusion would apply to t\ and % as shown in Table 2 



Table 2: Range of variation of parameters for S3 = — V3, q=9 



n 


n h 




ti (in units of tp) 


1 

(lo\) 4 (in units of Mp) 


tb (in units of tp) 


25 


16.8 to 19.2 


1Q -25 


6(10 16 ) 


10~ 9 to 4(10~ 9 ) 


h to 10 19 


20 


13.4 to 22.6 


10 -20 


2(10 13 ) 


10 -9 to 2(10~ 7 ) 


ti to 4(10 22 ) 


15 


10.1 to 25.9 


10~ 15 


10 10 


10~ 9 to 9(10~ 6 ) 


h to 8(10 25 ) 


10 


6.7 to 29.3 


10- 10 


5(10 6 ) 


10- 9 to4(10~ 4 ) 


h to 2(10 29 ) 


5 


3.3 to 32.6 


10~ 5 


2(10 3 ) 


10~ 9 to 2(10~ 2 ) 


h to 4(10 32 ) 


1 


0.7 to 35.5 


10- 1 


5 


10~ 9 to 0.4 


h to 2(10 35 ) 



3.2 Final results 

The results of the previous tables correspond to wavelengths that start inside the horizon and 
remain there, but there are other options. Depending on the value of S3 we can distinguish 
between 

(a) — 2 < S3 < — 1, in this region the physical wavelength decreases with time as long as 

t < t b . 

(b) — 1 < S3 < i, in this region the physical wavelength increases with time although not as 

fast as F _1 (t) 

In these two regions it is possible to distinguish between the following two behaviors: 

• the wavelength is outside the horizon at t = ti (initial time), which will occur for 
comoving wavelengths 



/ 4 \ 3 2 -s3 £3nl 



Then, the perturbation crosses inside the horizon at some time tnc- Although tuc can 
be both bigger of smaller than tb, only the case when is smaller will be of interest to us 
since we have been making the assumption t <^ t^. For a comoving wavelength Ao, the 
time of horizon crossing tnc is 
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tHC = 



Ao 



3o; \3 Jl + ssl i^a 
4 J 3 h 



The constraint tuc < % puts an upper bound on the comoving wavelength 



An < A n = 



|l-s 3 | \^oJ 6 



• the wavelength is inside the horizon at t = ti, that is Ao < Aj^. The perturbation will 
never cross outside the horizon since the slope for X(t) is negative at small times and 

2 i 

goes like fa at later times, while the slope of H grows always like t. 
3.3 Comparison with FRW Universe 

In this section we will compare the results that we have just obtained with those of an ordinary 
matte dominated FRW universe. The perturbations both inside and outside the horizon grow 
as [29], [30] 



5i - 



2 

t\3 



while the energy density changes as 



In analogy with the parametrization introduced in (3.1), let's parametrize the initial pertur- 
bations as 



8 = 10" 



2 

t \ 3 



and choose ti = Wtp, in the event that the universe starts out homogeneous and isotropic. 

3n 

The time at which the perturbations become of order one is t\ = 10~tp. 
Summarizing 



t FRW = 10 f tp (38) 

if 1 = lQi t P , (3.9) 



where 




at large wavelenghts 
at small wavelengths. 
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For the growth of perturbations for the Bianchi I model to be faster than in the FRW case: 
tf 1 < tf RW ; this constraint implies that the quantity s > |, which translates into the 
following bound for S3. 

S3 < 1 : for large wavelenghts 
S3 < — \ : for small wavelengths. 

Table 3: Comparison between the time that it takes for perturbations to grow to be of order one in 
the Bianchi type I model and in the FRW model, in the special case where S3 = — V3. The quantity 
n measures the size of the initial perturbations, n = 3 — log 5i for the FRW case, and n = s — logSi 
in the Bianchi case. 



n 


tf 1 (in units of tp) 


tf RW (in units of t P ) 


25 


6(10 16 ) 


3(10 37 ) 


20 


2(10 13 ) 


10 30 


15 


10 10 


3(10 22 ) 


10 


5(10 6 ) 


10 15 


6 


10 4 


10 9 


4 


5(10 2 ) 


10 6 


2 


20 


10 3 



4. Conclusions 

The Bianchi type I model (homogeneous and anisotropic) was introduced as a background 
spacetime for a cosmological model characterized by small density inhomogeneities. The evo- 
lution of these perturbations was investigated for a pressureless cosmic fluid and the equations 
turned out to be complicated by the presence of a coupling between the density contrast and 
the tensor perturbation modes. This coupling is absent for special choices of the anisotropy 
parameters and for isotropic models in general. Analytic solutions in the limits of small and 
large wavelength perturbations (compared to the horizon size) were found. The tensor mode 
that is coupled to the scalar energy mode turned out to be an oscillatory function of time 
with an amplitude decreasing in time for any possible value of the anisotropy parameter, 
whereas the density contrast showed a power law growing behaviour in time for a large range 
of parameter space. The growth is faster in the contracting direction of space compared to the 
directions that are in expansion for both small and large wavelength limits. These suggestive 
results motivated us to further analyze the growth of density inhomogeneities, specially, in 
view of the possibility that inflation might not occur if the spacetime is not homogenous. For 
quite a large range of parameters of the theory, the inhomogeneities in the energy density will 
grow to be of order one at a time where A is still a subdominant component of the energy 
density. As the period available for the fluctuations to become non-perturbative grows, while 
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the cosmological constant remains subdominant, the restrictions on the parameter space will 
weaken substantially. 

Since our analysis is restricted to perturbations, we cannot follow the evolution of in- 
homogeneities to the non-perturbative domain. In order to properly study the subsequent 
evolution of the universe, non-perturbative tools will have to be used. The results presented 
in this paper raise serious concern to motivate further work to find the fate of the universe 
in this context. 
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6. Appendix A 

The Kasner metric is generally written as 

ds 2 = -dt 2 + t 2pi dx\ + i 2pa dx\ + t 2p3 dxl (6.1) 
The coefficients pi, P2, P3 are constants satisfying the relation 

Pi +P2 +P3 =p\ +pl +pI = 1 (6.2) 

The three parameters can be equal in pairs in the cases (— 1,|,|) and (0,0,1). In all the 
other cases they are distinct, one being negative and the other two being positive. Under the 
assumption P3 < P2 < Pi, the allowed ranges follow from (6.2) 

2 

3<Pi<l, 

< P2 < §, 
~ < P3 < 0. 

The coefficients can be parametrized as 



Pi(u) ~ 

P2{U) ~ 



ps(«) - — 



l+u+u 2 ' 

1+M 
1+U+U 2 ' 

— u 
l+u+u 2 ' 



where u > 1 



The universe described by the metric (6.1) is spatially flat and, for any value of the co- 
efficients p,i, in expansion since the volume element is \J g( 3 ) d ? 'x = td 3 x. 

Prom (6.2) it is evident that the arguments pi of the scale factors in the different direc- 
tions of space cannot have the same sign, which indicates that the universe is expanding in 
some directions and contracting in some other directions with a law like = L^t Pi , where 
L? and L\ are respectively proper and comoving distances. 

A metric of the kind (6.1) can be generalized introducing a time dependence for the 
coefficients and can be put in the form 



ds 2 



--dr 2 + e 2a {e 2/3 ) ij dx i dx j 
7 



(6.3) 



where a = a(r), 0ij = (3(r)ij is a traceless diagonal 3x3 matrix. The function 7 = 7(2) can 
be eliminated redefining the time as dt = j~zdr. 

The equations for this metric, in the presence of an isotropic perfect fluid, were derived 
for the first time by [21] and later also by [24]. 



3 -1 

--dry - 37a - 37a 2 - jfyjPij = -{u + 3p) 

^dj + 7a + 37(d) 2 = ^(u-p) 
jfiij + 2jPij + fryafiij = 0, 

where the dot indicates the derivative with respect to r. It is easy to check that 



(6.4) 

(6.5) 
(6.6) 



*7 



b{jU -\- Cij 



solves the field equations if u acts as a new time coordinate 

„—3a 



du 



-dr 



(6.7) 



(6.8) 



with bij and Cy traceless diagonal constant constant matrices. The relation (6.2) can be 
translated into a condition involving the bij. The first of the conditions (6.2) translates into 
the requirement that Tr(/3jj) = 0, while the second condition (6.2) corresponds to 



Tr 



Sa + 2 



dfijj 
da 



d(3jj 
da 



9 



which is equivalent to 



El dfiii 

i=l,2,3 



6. 



(6.9) 
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Plugging (6.7) and (6.8) into Einstein equations, we have 



/ da 



1 



du 

Integrating, for p = au, we have 

_3(<r-l)a/2 = _ 



da , 



(6.10) 
(6.11) 



3^ Smh 



-b(l-a)u 



where b 2 = ^hjbij and ujq is an integration constant. For the special case of pressureless 
matter (a = 0), in terms of the proper time t 



„3a 



3u \ ( , b 
— - )t t + 4— 
4 / V w„ 



(6.12) 



7. Appendix B 

The perturbed metric is ds 2 = -±dr 2 + e 2a e 2 ^ (s!f + /ij(r, xfj dx l dx j . 
The unperturbed affine connection is 



TO _ 1 J30 dgoo 
1 00 — 2» dt 



00 — 2^ dt 
■p/c 1 „fcf rfffi 



0/ ~ 2# dt 





1 oo 



1 fci 



o 



a 9u r 



ki 



V _ l n 00c(9fei 
1 ki - 1,9 dt 



The first order perturbation in the affine connection is 



01 oo 



o ar& = o 

ST k. = Igkj &^ + l Sgkj dg» ^ = Q 



0i fci- 2^ at 01 ki - 2V y dx k a? a^ J 

The perturbed first order Einstein equations are 

«V = - ±6 9lu ,T - \g^{5g° x T aX + g aX 6T aX ). 
For a perfect cosmic fluid in a comoving frame we have [21] 



(7.1) 



- 7;jh - j(h + 2ah + 2ft, = <5u; + 3cSp 
2 ^ 
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djhi - d t k - 2Pitdjhb + 2{3 3t d 3 h\ + fcjdjh - p'abdiht = -^-{o J + P )5u l 

\e- 2a {e- 2 % s {d ia h a s + d sa h<* - {e 2 P) sk (e- 2 P) ab d ab h k 3 - d sj h) 
+ \ih) + )ph) + ^h(aSi + Pij) + + 7/^ - iPkA = \{8u - 5p)8) 

(7.2) 

where h = £ i=1>2)3 4 

Because spatial translation invariance is still a symmetry of the Bianchi I spaces, the 
modes can be separated by Fourier components. Let's select a particular wave-vector k, then 



h)(r, x) 


= ^{ry k - s 


5uj 


= W{T)e il * 


5p 


= P(r)e^ 


5u l 


= iVi(T)e il * 



(7.3) 



and /j, = Yli=i 2 3 A*i- ^ is useful to introduce a projection operator 

«ij = <>ij - tttt 

It S3jtlsfl6S . Ki^gKigj — ^ij — ^jii ^is &s — £111 d AC^s^s — •> W^llGrG / * /c — 0. 

We define the following parts of //'■ 



3 



s=l 

1 / t k s kt 



li = 

??} = /" S * ^j S /% - ^K ts Kij^ . (7.7) 

For simplycity we will consider a single component for the wave vector, k = (0,0,^3). 
After substituting (7.4)-(7.7) in (7.2) and doing some algebra, the equations (2.12) and (2.13) 
can be easily derived provided we make the identification: 

c W(T) ! 

= — --— and rj = r\ \ 

w(t) 
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where w(t) = Too- 

In the long wavelength limit with no residual isotropy (53 / ±2), the equations (2.12) 
and (2.13) can be combined into a fourth-order differential equation for 5, which can be solved 
to give 

a)) 

Substituting this result in the equation for r/, we have {t <C 

d 2 r] ldrj 2 -2(1+33) -2(i-"3) 1. ,ld<5 , 

For different choices of the parameters Sj and of the coefficients in (2.15), the behaviour 
of i] will be either increasing or decreasing as a function of time, but never oscillatory 

V = Ci + ^ + C 3 t 2J2 ^ (1 + C 4 \nt) (7.9) 

The constants Ci are related to the constants that appear in (7.8). 

The short wavelength limit translates into the condition FK 2 t 2 » 1. The trick to solve 
the equations in this limit is to average both 5 and r\ over an interval of time that needs to 
be relatively short in order to have a good approximation and, at the same time, long enough 
as to include several wavelengths. 

The behaviour in the short wavelength approximation in an isotropic model (where the tensor 
and the scalar modes are decoupled), is oscillatory with a decreasing amplitude (see Lifshitz 
and Khalatnikov, 1963, page 513). In the anisotropic case the situation is different because 
tensor and scalar modes are coupled, nevertheless an oscillating behaviour is expected in the 
very short wavelength limit, with a matter density perturbation smoother than the tensor 
mode. An average over several wavelength for r] would thus be negligible compared to an 
average for 5 if the former is represented as an oscillating function and the latter as a smoothly 
varying function. 

With this is mind, dividing both sides of the equation (2.12) by FK 2 t 2 and averaging, for 
t <^tb it simplifies as follows 

d 2 5 2(1 + s 3 ) dS 

~df + 3~t dt 
where now 5 is an average quantity. 
The solution to this equation is 

S = Ats- a -+B. (7.11) 

Notice that the density contrast has, on top of an overall oscillating behaviour which has 
been subtracted thanks to the averaging procedure, a power law changing amplitude which 
is positive for a large part of the range of possible values of the parameter S3. 
Replacing this into the equation for rj, we have 
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(7.10) 



cfrj 1 drj 2 



si - s 2 



A't 



-0.51 



(7.12) 



where the constant A' turns out to be of the order of 10 60 in units of sec -1 ' 49 according to our 
estimate of the parameters involved in perturbation theory, and we have assumed , S3 = — \/3. 
With a change of variable, we have 



X) + - + 10 4 x u 

X 



-Jl3.-0.51 



>' )A9 V 10 

where the dots indicate the derivative with respect to x = 



(7.13) 



j^, tp being the Plank time. 
The solution is a linear combination of Gamma, Hypergeometric and Bessel functions. In a 
simplified version of the previous equation setting the right hand side to zero (n is always 
bigger than one, so this is a reasonable approximation), the solution becomes 



77 = Ci J 



■ i 2 * 
hi -x4 



+ C 2 Y 



, 2 5 



(7.14) 



with C\ and C2 constants of integration and Jo and Yq are Bessel functions. This solution 
represents a sound wave with an amplitude that attenuates in time. 

Since the rigth hand side of equation (7.12) can be neglected without significantly altering 
the solution to the equation for r], we conclude that the choice of the anisotropy parameters 
Si doesn't really affect the overall behaviour of the tensor mode. This is radically different 
from the scalar mode, which is dramatically affected by that choice (7.11). 
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